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A toy model for quantum spin Hall effect 
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In this Communication, we investigate a toy model of three-dimensional topological insulator 
surface, coupled homogeneously to a fictitious pseudo spin-| particle. We show that this toy model 
captures the interesting features of topological insulator surface states, which include topological 
quantum phase transition and quantum spin hall effect. We further incorporate an out-of-plane 
magnetic field and obtain the Landau levels. 


PACS numbers: 73.43.Nq,73.43.-f, 85.75.-d,75.45.+j 

Introduction.- In recent years, topological insula¬ 
tors (TIs) have captivated considerable attention of re¬ 
searchers [1-5, 8]. These fascinating materials involve 
compounds such as Bi 2 Se 3 or Bi 2 Te 3 , whose electronic 
bulk structure is an insulator with a finite gap separat¬ 
ing the conduction band and the valence band, but their 
edges (for 2D TIs) or surfaces (for 3D TIs) have gap¬ 
less states, which are protected by time-reversal sym¬ 
metry (TRS). These states are robust to perturba¬ 
tions that do not break this symmetry. In recent years, 
these gapless states have been the focus of interest due 
to the simple Dirac-like Hamiltonian which describes 
them. However, breaking TRS introduces many inter¬ 
esting phenomena; this can be achieved by depositing 
a ferromagnet or a superconductor on the surface of a 
TI. Thus, the topologically protected surface state de¬ 
velops a gap, which leads to interesting electronic trans¬ 
port, which include half quantized conductivity, Majo- 
rana bound states, etc.fl, 11-15, 17, 18]. They also 
have potential technological applications in the field of 
spintronics[19, 20]. 

In this Communication, we study a toy model that cap¬ 
tures some of the interesting features of surface states 
in topological insulators. Specifically, we study a ficti¬ 
tious pseudo spin-1 particle interacting homogeneously 
with the surface of a TI. Assuming anisotropic “in-plane” 
and “out-of-plane” exchange interactions on the interface, 
say A|| and A^, the coupled system exhibits a topologi¬ 
cal quantum phase transition as a function of ^ = X±/X\\. 
This phase transition is associated with a quantum phase 
transition point at ^ = 1, which separates two regions; 

with two topologically protected Dirac points 
(semimetallic phase) and ^ > 1, which is fully gapped 
(quantized spin Hall phase). We further show that when 
the ^-component of the pseudospin is conserved, the sys¬ 
tem decouples into two Hamiltonians which are related 
by TRS. The coupling term opens two gaps at the T 
point on the interface, which are degenerate but with op¬ 
posite spin orientations. Each gap opening gives rise to 
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a half quantized Hall conductivity. Upon applying an 
electric field, an opposite spin current is induced on each 
interface leading to a quantized spin Hall conductivity. 
Furthermore, we introduce a non-conserving term in the 
quantized spin Hall Hamiltonian by breaking the conser¬ 
vation of the ^-component of the pseudospins. This leads 
to the obliteration of the quantized spin Hall conductiv¬ 
ity. However, in this non-conserving regime we obtain a 
nontrivial topological spin Chern number by diagonaliz¬ 
ing the projection of the pseudospin onto the occupied 
valence bands [21]. 

Simplest Single Dirac Point- Most of the interesting 
physics of 2D metals are captured by investigating the 
robustness of the band touching points (nodes). The 
bulk energy band near these points usually replicates a 
2D Dirac Hamiltonian; thus these nodes are called Dirac 
points. The simplest form of Dirac point can be mapped 
out from this simple Hamiltonian: 

H = vf{z X ct) • k, (1) 

where vp is the Fermi velocity, ai,i = x,y, z are the Pauli 
matrices representing the real spins of the surface states. 
The ^-component is taken to be perpendicular to the 
plane of the 2D sample, and k = ky) is the 2D surface 
Brillouin zone momentum. This Hamiltonian describes 
many known physical systems, such as the surface of a 3D 
topological insulator [8] , and the Dirac point Hamiltonian 
in graphene in which the Pauli matrices are pseudospins 
[10]. Evidently, this Hamiltonian is invariant when the 
spin and momentum go to minus of their original values; 
in other words, the Hamiltonian possesses TRS, where 
the time reversal operator is given 0 = iayJC; /C is a 
complex conjugation. The energy spectrum is trivial and 

given by Ss = svF^k‘^ + ky', s = T. It is apparent that 

the two bands touch each other at k = 0. The robustness 
of this Dirac point is guaranteed by any perturbation 
that preserves the TRS of Eq.(l). By breaking TRS, 
the degeneracy of the band might be lifted and many 
interesting phenomena can emerged. There are several 
terms that break TRS; for instance, Axax and Ay ay do 
break TRS; they correspond to Zeeman magnetic field 
contributions along the x and y directions respectively. 
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However, addition of these terms to Eq.(l) simply shifts 
the position of the Dirac points to k = (0, —Ax/vf) and 
k = {Ay/vp, 0) respectively. Thus, the degeneracy is not 
lifted. The only TRS breaking term that lifts the degen¬ 
eracy of the bands is the direct coupling of the surface 
electron spins to the Zeeman magnetic field term per¬ 
pendicular to the plane, i.e., A^a^. This term can also 
be generated by depositing a ferromagnet on the surface 
electrons [11, 13, 14]. It is evident that this contribu¬ 
tion opens a gap of size 2\Az\ at k = 0. Consequently, 
when the Fermi energy lies between the gap, this leads 
to a half-quantized Hall conductivity [8, 9]. In graphene, 
however, there are an even number of Dirac points. In 
this case, the charge Hall conductivity vanishes in the 
ordinary insulating state [5-7]. 

Toy Model- As mentioned in the preceding sections, 
the surface of a 3D topological insulator possesses many 
interesting features when TRS is explicitly broken. In 
this section, we will consider a 3D topological insulator 
whose surface is homogeneously coupled to a fictitious 
pseudospin-^ particle; the low-energy effective Hamilto¬ 
nian can be written as 

TL = vf{z X tr) • k - \\\Tx(7x - ^^TzCFz, (2) 

where (Ji and r^, {i = x^y^z) are the Pauli matrices 
acting on the topological insulator space and the ficti¬ 
tious pseudospin-1 particle space respectively. The cou¬ 
pling constants Ay and are the anisotropic in-plane 
and out-of-plane homogeneous exchange interactions. In 
most cases of physical interest, such interaction is usually 
inhomogeneous, z.e., the last two terms in Eq.(2) should 
contain a delta function. Thus, Eq.(2) does not describe 
any known physical system. However, it is possible that 
it might be applicable to a pseudo quantum qubit, but 
we are not interested in any specific system, instead we 
will consider it as a toy model that captures some of the 
interesting physics of TI surface states. Although our 
model does not describe any known physical system, it 
possesses some features that are similar to other models 
that describe known systems. This is the main purpose of 
this Communication. We will consider Eq.(2) as the basis 
of our investigation in this Communication. It is appar¬ 
ent that Eq.(2) explicitly breaks TRS^ i.e , [0,H] ^ 0, 
where O = Tx ^icFylC. Due to the fictitious pseudospin-^ 
particle, Eq.(2) is obviously a 4-band model. 

Topological Quantum Phase Transition - Most known 
systems such as graphene and thin film topological in¬ 
sulators possess topological quantum phase transition. 
This is manifested at the gap closing point which sep¬ 
arates an ordinary insulator from a quantized Hall insu¬ 
lator [16, 26]. These phases are usually distinguished by 
a topological invariant quantity called the Chern number. 
It is generally defined as [25] 

JSkgdk^g-^Qdk,g-^gdk^g-\ (3) 

where is the totally antisymmetric tensor, a = 

0, X, etc., labels the components of a four-vector, and 


the Matsubara Green’s function can be written as 

g~^{iuJn, k) = iWrJ-ixA “ "H. (4) 

The identity I 4 x 4 is a 4 x 4 matrix and = {ioon, k) 
is the momentum four-vector. In principle the Chern 
number can be computed for any model of interest. One 
finds that it has a unique value in each phase, which 
is immutable by any smooth deformation of the system 
(provided the gap does not close). In general, the emer¬ 
gence of quantum phase transition in electron systems 
requires the violation of CPT symmetry [25], i.e., charge 
conjugation, parity, and time-reversal symmetries. It is 
evident that our model in Eq.(2) explicitly breaks parity 
and time-reversal symmetries. Thus it can capture quan¬ 
tum phase transitions similar to those predicted in known 
systems. In order to see this, we diagonalize Eq.(2) and 
find that the eigenvalues are given by 



where s = A and r] = 0,1. There are two energy sectors 
Spy and £-y. In each sector there are two bands with 
7^ = 0 being the conduction band and r] = 1 being the 
valence band. The £-y bands have two Dirac points 



FIG. I. Color online. The phase transition in space 

with = I = A||. The quantum phase transition point ^ = I 
separates two regions; ^ < I, with two Dirac points and ^ > I, 
which is fully gapped; see Fig. (2). 

bands are always gapped) at ±K, where 

K= (o, ^ = Ax/A||. (6) 

The topological quantum phase transition as a function 
of ^ can be understood as follows. At the point ^ = 0, the 
system is a semimetal with two Dirac points separated 
by the wave vector vpky = 2Ay along the ^-direction. 
The semimetallic phase is controlled by the in-plane ex¬ 
change interaction. Provided ^ < I, these two Dirac 
points remain intact, separated in momentum space by 
a wave vector vpky = 2A||y/l — They are topolog¬ 
ically protected and cannot be eliminated by changing 
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the parameters of the Hamiltonian. However, as ^ is 
increased, the two Dirac points remain stable and sub¬ 
sequently annihilate each other at the phase transition 
point ^ = 1, and emerge as a single Dirac point at the F 
point K = (0,0). The gap further reopens as one moves 
to the region ^ > 1 and the system becomes an insulator 
with a gap of 2A|| (^—1) at K = (0, 0). As shown in Fig.(l) 
and Fig. (2), the quantum phase transition point ^ = 1 
separates two regions; ^ < 1, with two Dirac points and 
^ > 1, which is fully gapped. The homogeneous in-plane 
coupling AII plays a similar role as the parallel magnetic 
field in a thin film topological insulator [26]. 



FIG. 2. Color online. The energy bands along the ky direction 
for different ranges of the parameters; see text for explanation. 


£ 



FIG. 3. Color online. Schematic diagram of the two Kramers 
energy bands on the interface with a gap of 2A at the F point. 

where Xs(A) = ^ (^^1 + s-^, -ise^^’= - s^) and 

Ok = tan“^ (ky/kx). Fig. (3) shows the Kramers pairs 
of the energy bands for spin-up and spin-down, with an 
energy gap of 2A at the F point {kx = O^ky = 0). At zero 
temperature, and when the Fermi energy lies between the 
gap, the conduction band is empty and the valence band 
is full; the Hall conductivity is given by [29] 


a 


Tz — 
xy 



( 11 ) 


where 


Quantum Spin Hall Effect- Quantum spin Hall ef¬ 
fect is known to be manifested in many known systems 
[2, 22, 27, 28, 31]. In this section we will show how the 
toy model in Eq.(2) captures this interesting physics. Let 
us consider the limit Ay —?► 0. In this limit there are no 
Dirac points, the system is fully gapped. The Hamilto¬ 
nian decouples into two copies which are related by TRS: 

«=«*)' p) 

where 

y-U = vrikycr^ - k^ay) ± Xa^, ( 8 ) 

and A = A^. The eigenvalues are pair degenerate (r^ is 
conserved), given by 




it 


— I d^kTr., 


( 12 ) 


is a topological invariant Chern number, = 

(V X is the Berry curvature and = 

i is the Berry connection. Expressing the 

differential operator in polar coordinates yields: 


1 ( d{kAe,) dAk\ 
k \ dk dOk ) ’ 


(13) 


where A 0 ^ = (■*/'!? |(90^V’-) and = i (V’l* )• 

Since Ok appears as a U{ 1 ) phase in the eigenspinors, 
Aqj^ and Ak are independent of Ok, thus the second term 
in Eq.(13) vanishes and the resulting expression becomes 
[30]: 


Sg = sfvfH-Afi = sek, 


(9) 


where k = yk^Eky. The corresponding eigenspinors 
are found to be 



= M^^ik) = i{rj\d0,rA}, ( 14 ) 

Integrating Eq.(12) using Eq.(14) yields a half quan¬ 
tized Hall conductivity at each gap opening: 

-e^ 

= ('r-5)^sgn(A). 


( 15 ) 
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Evidently, the Hall conductivity is coupled to the ficti¬ 
tious pseudospin orientations; see Fig. (3). Thus, an ap¬ 
plication of an electric field will induce an opposite ficti¬ 
tious spin current. The charge-Hall conductivity vanishes 
due to TRS: 


= 0 , (16) 

while the resulting half quantized spin Hall conductivity 
is nonzero: 




- Ye 




It should be noted that this spin Hall conductivity is not 
measurable. Physical known systems in which spin Hall 
effect has been realized include HgTe/CdTe semiconduc¬ 
tor quantum wells [ 22 ], and the ultrathin film topological 
insulator [27]. Thus our coupling constant A plays a sim¬ 
ilar role as the coupling constant in these systems. 

Effects of Non-Conserving Terms and Magnetic Field- 
Our analysis for the quantum spin Hall effect in the pre¬ 
vious section relied on the fact that the ^-component of 
the fictitious spin 1/2 system is conserved. We will now 
investigate the effect of adding a non-conserving term to 
the Hamiltonian. There are many ways to do this, but 
we will choose the most plausible way. We assume that 
the fictitious particle is a two-level system with degen¬ 
erate minima. Tunneling obviously breaks this degener¬ 
acy, which introduces dynamics into the system and thus 
breaks the conservation of the ^(-component of the spin. 
This dynamics effectively corresponds to a term of the 
form AxTx. Thus, the model in Eq.(7) then becomes: 


n = 




(18) 


In this Hamiltonian, the off-diagonal term is akin to 
structure inversion asymmetry potential in a thin film 
topological insulator [31]. Diagonalizing this Hamilto¬ 
nian we find that the eigenvalues are given by 


Ssr, = {-irf^Fk + sA.f + X^. (19) 

The corresponding eigenspinors for the conduction 
band (77 = 0 ) and the valence band (77 = 1 ) are given 
respectively by 

( 20 ) 


Since is no longer conserved, the Hall conductiv¬ 
ity in Eq.(15) is obliterated. However, in the regime 
^ I 5 0116 can still define a pseudo spin Chern num¬ 
ber by diagonalizing the operator VtzV [21, 30], where 
V = XsXl is the projection operator onto the occu¬ 

pied valence bands in Eq.(21). The matrix elements are 
expressed as {Xs\^z ^ Icr|Xs ')5 s^s' = ±, and the eigenval¬ 
ues of this matrix are given by uj±{k) = ±|p(/c)|, where 

p(/c) = 2 ( 7 X++ 7 X_ — u _ U-\ _), (23) 

and 77 = 1 always. The corresponding eigenfunctions of 
Tz are given by 

=/|(x+±x-)- (24) 

The pseudo spin Chern number is now given by Eqs.( 12 ) 
and (14) by replacing ipR with in A4{k). Computing 
the corresponding inner products yields: 

M±{k) = ±^ + l. (25) 

Integrating Eq.( 12 ) one obtains the Chern number: 

= Al*(oo) - A4*(0) = (p(oo) - p{0)), (26) 

which defines a nontrivial topological spin Hall phase pro¬ 
vided that Ax/X 1. The function p{k) is depicted in 
Fig. (4) for a specific value of Ax/X. It is analogous to 
the Pfaffian [2, 31], which is obtained by diagonalizing 
the TRS operator 0 in the occupied valence bands ba¬ 
sis. 



( 21 ) 

where 

fT,- 


FIG. 4. Color online. Plot of pfk) vs. k with Ax/X = 0.04. 

Interesting features also arise in the presence of a mag¬ 
netic field. Now let us apply a constant magnetic field 
along the out-of-plane direction in the system described 
by Eq.(18). The resulting Hamiltonian becomes: 

H = vf{z X cr) • kA + (A^ - - A^t^, 


1 

Usr]± — - 


(27) 
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where Ica = (—iV + ^A) is the magnetic field experi¬ 
enced by the surface electrons through the Peierls sub¬ 
stitution, also known as the orbital effects on the elec¬ 
trons. Here A = —yBx is the vector potential for the 
Landau gauge. Diagonalization of Eq.(27) is achieved by 
introducing the creation and annihilation operators: 

6 = -%(n, + zn,); = (28) 


where = cjeB is the magnetic length and 11 = —iV + 
^A. Thus, Eq.(27) becomes: 

1-L = iujBV^icr^b - + (A;^ - - A^r^, (29) 

with ujb = vfIIb being the characteristic frequency anal¬ 
ogous to cyclotron frequency. The corresponding landau 
levels for n 7 ^ 0 are obtained as: 

f^± = 2c^|n + A^ + A^,-FA^ 

±2^2uj%nAl + Al (A^+A^), (30) 


There are two energy sectors for each Landau level which 
are controlled by several parameters. Both sectors be¬ 
have quite differently as a function of the magnetic field. 
The Hall conductivity can as well be computed using 
these Landau levels [23, 24]. 

Conclusion- In conclusion, we have studied a toy 
model that captures the interesting physics of quantum 
spin Hall effect and topological quantum phase transition 
on the surface of a 3D topological insulator. Our model is 
comprised of a homogeneous exchange coupling between 
the surface electrons and the fictitious spin-^ particle; 
we showed that this model exhibits both quantum phase 
transition and half quantized spin Hall effect. We pre¬ 
sented a lucid exposition of the quantum phase transition 
and the topological invariant quantity that characterizes 
the Dirac points. Introducing dynamics in the fictitious 
particle leads to a non-conserving term in the Hamilto¬ 
nian, which breaks the conservation of ^-component; thus 
obliterates the spin Hall effect. We obtained a nontriv¬ 
ial topological spin Chern number by diagonalizing the 
projection of the pseudo spin onto the occupied valence 
bands. Eurthermore, we obtained the Landau levels in 
the presence of a magnetic field. 
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